The Topological Non-connectivity Threshold and magnetic phase transitions in 
classical anisotropic long-range interacting spin systems 
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We analyze from the dynamical point of view the classical characteristics of the Topological 
Non-connectivity Threshold (TNT), recently introduced in F.Borgonovi, G.L.Celardo, M.Maianti, 
E.Pedersoli, J. Stat. Phys., 116, 516 (2004). This shows interesting connections among Topology, 
Dynamics, and Thermo-Statistics of ferro/paramagnetic phase transition in classical spin systems, 
due to the combined effect of anisotropy and long-range interections. 
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I. INTRODUCTION 

The magnetic properties of materials are usually de- 
scribed in the frame of system models, such as Heisenberg 
or Ising models where rigorous results, or suitable mean 
field approximations are available in the thermodynami- 
cal limit. On the other side, modern applications require 
to deal with nano-sized magnetic materials, whose intrin- 
sic features lead, from one side to the emergence of quan- 
tum phenomenaT], and to the other to the question of 
applicability of statistical mechanics. Indeed, few parti- 
cle systems do not usually fit in the class of systems where 
the powerful tools of statistical mechanics can be applied 
at glance. In particular, an exhaustive theory able to 
fill the gap between the description of 2 and 10^^ inter- 
acting particles is still missing. Moreover, also impor- 
tant well-established thermodynamical concepts as the 
temperature, become questionable at the nano-scale0. 
In a similar way, long-range interacting systems belong, 
since long, to the class where standard statistical me- 
chanics cannot be applied tout court. Indeed, they dis- 
play a number of bizarre behaviors, to quote but a few, 
negative specific heat jSj and hence ensemble inequiva- 
lence 0, temperature jumps, and long-time relaxation 
(quasi-stationary states) l^. Therefore, from this point 
of view, few-body short-range interacting systems share 
some similarities with many-body long-range ones. 

Within such a scenario, and thanks to the modern com- 
puter capabilities, it is quite natural take a different point 
of view, starting investigations directly from the dynam- 
ics, classical and quantum as well |l2j| . It was in this spirit 
that quite recently in a class of anisotropic Hciscnberg- 
like spin lattice systems, a topological non-connection of 
the phase space was discovered JlSj. Initially, for histor- 
ical reasons @ , this was referred to as broken ergodicity, 
since if the phase space is decomposable in two uncon- 
nected parts (hence topologically non-connected) then a 
breaking of the ergodicity is indeed a trivial consequence 
0. Here we prefer to call Topological Non-connectivity 
Threshold (TNT) the value Etnt where such a disconnec- 
tion sets in upon lowering the total energy E of the sys- 
tem. This result was found, first numerically and later 



analytically, in a class of spin models where important 
and rigorous results have been obtained during the last 
century, though generally only in the thermodynamical 
limit (total number of microscopic constituents of the 
system N ^ (x keeping fixed the number density N/W^, 
where R is the spatial size of the system and d is its spa- 
tial embedding dimension). Nevertheless, to the best of 
our knowledge, apparently nobody took care of the dy- 
namics, and consequently nobody spoke of this simple 
but relevant property. 

This dynamical point of view has a few interesting 
classical consequences. First of all it explains, from the 
point of view of microscopic dynamics, the possibility of 
ferromagnetic behavior in small system. Indeed, in ab- 
sence of external field and external noise (temperature) a 
magnetized system, (belonging to one branch of the non- 
connected phase space) remains magnetized simply be- 
cause it cannot move to the other one. Furthermore, our 
TNT is surely related to recent results |S 1^ ITol| connect- 
ing topological transitions (TT) and thermo-statistical 
phase transitions (FT), even if such investigations again 
concern the thermodynamical limit only, and they re- 
late to usual FT of canonical thermo-statistics. How- 
ever, it has been recently stressed 111] that microcanon- 
ical thermo-statistics is the theoretically more suitable 
description for systems with small size and/or long-range 
interactions. 

We therefore begin in Section ^1 with a short descrip- 
tion of our class of models and the topological properties 
of the TNT for finite and infinite N, pointing out the cru- 
cial role played by the first key ingredient of our models, 
the XY anisotropy In particular we would like to answer 
the following question: what happen to the energy Etnt 
and to the corresponding iV-spin configuration when 
the N oo? It is exactly answering this question that 
the deep connection between TNT and the other key in- 
gredient, the long-range interactions, becomes apparent. 
Our results 0| can be summarized as follows: 

1) For finite sized systems TNT exists for short range 
and long range as well, as soon as anisotropy persists. 

2) When the number of particles becomes large the 
energy size of the non-connected region goes to zero with 
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respect to the total energy size for short range interaction 
while it goes to some finite quantity for long range ones. 

In Section IIIII we switch to the dynamical properties 
found in 15] which accompany such a special (and rather 
"big") TT, and its relations with its thermostatistical 
properties, namely the occurence and location of a usual 
(paramagnetic/ferromagnetic) PT, using techniques from 
large deviation theory 16] within the microcanonical de- 
scription of the system. 

We conclude in IIVI Everywhere here we restrict our- 
selves to the classical case; for a recent discussion of its 
quantum counterpart, we refer to p^llSj. 



effects are conveniently encapsulated within two simple 
and easy-to-handle parameters a and 77, in order to make 
the study of the minimum-energy configurations analyt- 
ically tractable at iV — > 00 and numerically feasible at 
finite N. Concerning the physical motivations for such 
choices, we refer to the recent discussion in [Tsj . 

Depending on the specific A'^-spin configuration the 
corresponding energy E = H{S^) will vary according 
to 10. One (not necessarily unique) configuration S^^^, 
to be specified soon, turns out to have minimum energy 
Emin, defined by the minimum of the Hamiltonian (|T|l 
over all conceivable configurations: 



II. MECHANICS AND TOPOLOGY 



Emin ^ Min[H] 



(3) 



As a paradigmatic model example of TNT, let us con- 
sider the following class of lattice spin models, described 
by the Heisenberg-like Hamiltonian: 



(1) 



where Sf,Sf,Sf are the spin components, assumed to 
vary continuously; i,j = 1, ...N label the spins positions 
on a suitable lattice of spatial dimension d, and is the 
inter-spin spatial separation therein. Here for simplicity 
we consider a d = 1 lattice. (See jlj] for extensions to 
d ~ 2,3.) For later notational convenience, we define 
the single-spin vector S = {Sf,Sf,Sf) and the A^-spin 

configuration — {81,82, ■■■Sn)- The tip of each i-th 
spin is taken to lie on the unit sphere, i.e. it satisfies 
the constraint (Sf)^ + {Sff + {Sf)'^ = 1. Also, a > 
parametrizes the range of interactions (decreasing range 
for increasing a) and < 77 < 1 parametrizes the XY 
anisotropy. For a = 00 we recover nearest-neighbor in- 
teractions, while a — corresponds to infinite-range in- 
teractions. A mean-field model is obtained by setting 
a = and including as well the (somewhat non-physical) 
self- interaction pairs i ~ j'. 
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H,nf^jN'imxy--N'imy) 



(2) 



where mx,y,^ = {^/N)J2, S^'^'^ While this might be 
thought of as a negligible modification for ^ 00, 
nevertheless it has non-negligible effects concerning the 
chaoticity properties of the system. Indeed, the dynam- 
ics of the truly mean-field a = system turns out to 
be exactly integrable'l5]. Here we are not interested in 
the most general spin Hamiltonian giving rise to a TNT 
(for instance in [isL lislj a term containing a transversal 
magnetic field Bz has been added to H). Rather, we 
focus on the very simple Hamiltonian which already 
contains the two essential ingredients which give rise to 
the TNT, i.e. anisotropy and long-range, whose main 



Since < 77 < 1 the minimum energy configuration S^j^ 
is attained when all spins are aligned along the Y axis, 
which defines implicitly the easy axis of magnetization. 
In the same way, let us define the TNT energy Emt as 
the minimum energy compatible with the constraint of 
zero magnetization along the easy axis of magnetization: 



Etnt = Min [H \my^O] 



(4) 



corresponding to some other A^-spin configuration 
to be specified later. By definition, in general Ej^in l£ 
Etnt, and in particular it may be that Emin 7^ Etnt- We 
call this situation is topological non-connection, as will 
become clear in a moment, and its physical (dynamical 
as well as statistical) consequences are rather interesting. 
Indeed, consider a system prepared at time t = with a 
definite sign of magnetization, say my > and an energy 
value Emin l£ E < Etnt- As time goes by, the system 
evolves upon the constant energy surface H{S'^) = E in 
configuration space. Nevertheless, due to the continuity 
of the dynamical equations of motion the magnetization 
my{t) (not a constant of motion) may well change its 
size, but it can never change its sign, instead. Indeed, 
in order to assume a value rriy < it should have to go 
through at least one configuration with rriy = 0, which 
by definition cannot belong to the E < Etnt surface. The 
whole situation can be summarized as follows. 

Topology: in configuration space the surface at fixed 
energy E is topologically non-connected in two com- 
ponents, each characterized by a magnetization either 
TTiy < or niy > 0. 

Dynamics: though the two components are energeti- 
cally accessible on equal grounds, the ergodicity of the 
constant E surface is trivially broken, since there exist 
no dynamically allowed path inbetween them. 

Thermo- Statistics: de-magnetization is in principle im- 
possible below the TNT, so we may speak in some sense 
of a ferromagnetic phase. Of course, the application of 
a magnetic field, or a thermal noise, can give the en- 
ergy necessary to overcome the energy barrier, thus in 
principle - if not in practice - allowing for a magnetic 
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reversal. On the contrary, for energy values E > Etnt, 
de-magnetization is in principle possible, and we may 
speak in some sense of a paramagnetic phase. However, 
being above the TNT docs not automatically guarantee 
that, for any combination of parameter values and initial 
conditions, a system initially magnetized at t = will 
for sure eventually de-magnetize within a given finite ob- 
servational time tabs- First, as reported in ,15.] some in- 
variant multidimensional structures can appear in phase- 
space preventing the motion from covering the whole en- 
ergy surface. This lack of ergodicity is well known in 
2D Hamiltonian systems, where KAM tori prevent mo- 
tion from wandering the whole phase space |l9j| . while 
it turns out to be more complicated in multidimensional 
system, see for instance the famous Fermi-Pasta-Ulam 
problem (for interesting reviews see [20l l2ll|'). Such in- 
variant structures usually disappear in the presence of 
strongly chaotic motions 19]. We can therefore say that 
strong chaos is somehow another necessary ingredient in 
order for the system to be in its paramagnetic phase. Sec- 
ond, even given strong enough chaoticity to "encourage" 
the system to explore all the available phase (or configu- 
ration) space, yet the system could be given not enough 
time to actually do it, so effectively "freezing" it within 
the niy component where it started from. 

For finite N systems the XY anisotropy is the only 
necessary ingredient in order to have Etnt < Emm and 
hence the TNT. For N ^ oo, one quickly realizes that 
Emin — > — oo and guesses that Etnt —^—oo, but may still 
wonder wehther Etnt — > ^-mm as well, thus making the 
TNT physically irrelevant. So we define the topological 
non-connection ratio: 



Etnt — En 

Eniin 



(5) 



which expresses how large a fraction of the energy range 
is topologically non-connected. Here we use the denomi- 
nator Emin (instead of the whole energy range), since, in 
such kind of models, one is usually interested in the neg- 
ative energy range only. Correspondingly, we will refer to 
a system as topologically non-connected if r ^ const ^ 
in the limit N —>■ oo. 

To show why and how this whole comes by in our 
models, we now concentrate on the energy E of the 
above mentioned interesting A^-spin configurations. The 
minimum-energy configuration (all spins aligned 

along the F-axis), can also be thought of as composed 
of 2 equal blocks of N/2 spins, all up. Correspond- 
ingly, the minimum energy Emin < can be split into 
3 contributions, namely 2 equal intra-block bulk energies 
i?l < plus 1 inter-block interaction energy V|| < 0. 
Flipping just 1 block, its bulk energy does not change, 
i.e. El = El < 0, but the interaction energy reverses sign 
to Viji = —Vii > 0. However, in such a new configuration 
S^l the magnetization is now iriy = 0. Both bulk energies 
are as low as they can, while the constraint rriy = frus- 
trates the total energy to rise above Emm by an amount 
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FIG. 1: TNT spin configuration. Here is a = 0.1, iV = 100, 
rj — 0.9. Only the central part of the chain has been shown. 



2|V||]. We therefore take E^ = Emin + '2\V^^\ as a some- 
what obvious candidate for Etnt and the related spin con- 
figuration S"!^ as the TNT spin configuration S^j. 
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1 example of a numerically obtained S^, is given in 



Fig.n As one can see, excluding a small central domain, 
the spins are arranged in two equal-size blocks with op- 
posite spins, i.e. the real S^t is remarkably similar to, 
though in fact somewhat slightly different from, the ideal 
S^^. Extensive numerical simulations using constrained 
optimization and analytical estimates as well 14], confirm 
such an impression, with natural analogues in d = 2,3. 

Quite generally, the exact TNT configuration at finite 
N for arbitrary a and 77 can be found only numerically. 
Nevertheless, concerning the large N limit it is possible 
to make a few definite statements First, if a ^ 0, 



rj — 1, and A^ sufficiently large, then S 



cN 



the 



twofold sense that the TNT spin configuration is approx- 
imately half up and half down, and the energy difference 
Eu due to the domain is a finite quan- 



E, 



Et, 



tity independent of A^. Second, for A^ — > 00, ii a > d 
(short-range) then r ^ but if a < 0? (long-range) then 
r — > const 7^ 0. In particular, if d = 1 then Too = 2 — 2". 
Intuitively, while for short-range the bulk volume contri- 
bution El + El eventually dominates over the geometric 
inte(sur)face area energy interaction V||0, for long-range 
interactions the "effective interface" is the whole bulk as 
well. As now both the inter-block and the intra-block in- 
teraction involve all spin pairs in the two blocks, though 
each pair with a different intensity oc energies will 
then scale like the number of pairs oc A^ times the in- 
teraction at typical distance Vij ~ R. So there holds 
the same scaling \Vii \ ^ {R^f/W - A^-"/"* - \Eminl 
though with different (and a-dependent) proportionality 
constants, and so r ~ const 7^ as A^ — > 00. 



III. DYNAMICS AND THERMO-STATISTICS 

Here, we show how to impose a self-consistent Hamil- 
tonian dynamical evolution upon the spin systems de- 
scribed by Afterward, we will employ such dynamical 
equations, numerically integrated for long enough times, 
to study the time-evolution of S^{t). Following 0,0] 
here we focus on the time evolution instantaneous mag- 
netization nriy(t) started with some E and ^^(O) and 
evolved according to jnjl. Complementarily, we look at 
its statistical distribution Pens {E, my), built via a ran- 
dom sampling of an ensemble of initial conditions, all 
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FIG. 2: Left column (a,b,c) : time evolution of the magnetiza- 
tion my{t) for different specific energy values. Right column 
(d,e,f) : probability distribution of the magnetization at a 
given specific energy. Parameters are a = 0, r; = f , A'^ = fO. 
Upper line ((a,d) below the TNT ) : E/N = -0.7. Middle 
line ((b,e) between the TNT and the statistical thresholds 
E/N = —0.3. Lower line ((c,f) above the statistical thresh- 
old) E/N = O.f. 



with the same E, tunable at will. Finally, we will show 
a connection between Dynamics and Thermo-Statistics. 

For each i the spin components Sf,Sf,S^ satisfy 
the usual commutation rules for angular momenta, i.e. 
{Sf, Sj} — SijSf, (and cychc) where {, } are the canon- 
ical Poisson brackets. As usual, (and this procedure im- 
mediately translates to the quantum case |17l ) , start- 
ing from the Hamiltonian (^) we straightforwardly derive 
the dynamical equations as: 



{H, S,} 



(6) 



As is well-known, for such dynamical equations the total 
energy E and the spin moduli 15^1^ = 1 are constants 
of the motion. Usually, for energies E not too close to 
Emim the dynamics is characterized by a positive largest 
Lyapunov exponent, corresponding to strongly chaotic 
motion [l^ . On the contrary, the dynamics of the truly 
(all pairs ij) mean-field a — system turns out to be 
exactly integrable'l5]. Given strong enough chaoticity, 
typical examples of evolution curves of my{t) and corre- 
sponding Pens{E,my), are shown in Fig.|5| 

Again, the whole situation can be summarized as fol- 
lows. 

Thermo-Statistics: for E < Emt (Fig. the two 
peaks of the distribution are well separated, while for 
E > Etnt (Fig. Et) they come close. At Estat they merge 
into one single peak (Fig.|2f), as expected from the sta- 
tistical analytical estimate. Note that this simple (and. 




(E-E„)/N 

FIG. 3: Divergence of magnetization reversal times close to 
the TNT. Here is a = 0, = 1, and different A'' values as 
indicated in the insert. Lines are the best fit according to 
Eq. Q. 



in its roots, purely geometrical) argument immediately 
prooves that, in general, Etnt < Estat- excluding catas- 
trophic situations, upon increasing E the two separated 
inner tails must first touch at Etnt before the relatively 
outer twin peaks may eventually merge at E^tat- Indeed, 
this very same argument automatically applies to any 
model with probability distribution changing from sin- 
gle to double-peaked , e.g. the discrete mean field ip^ 
model 0. In other words, if both a TNT and a FT are 
known to exist, then in general they will occur at differ- 
ent energies. And depending on wether E is increased or 
decreased, both a TNT and a FT in some sense "antici- 
pate" each other. Of course, for a given specific system 
nothing automatically guarantees that either a TNT or 
a FT do indeed both exist. So, while in principle a TNT 
and a FT neither imply nor exclude each other, in prac- 
tice they "strongly suggest" each other. 

Dynamics: at low energy (a) the system is always mag- 
netized with small fluctuations, and at high energy (b) 
the system is on average non-magnetized with large fluc- 
tuations. In between (c) the magnetization jumps errat- 
ically up and down. We can usefully define a magnetic- 
reversal time-scale |l5f . as the average time necessary 
to magnetization to reverse its sign. In the presence of 
strong chaos (dependent on the parameters N and E) 
magnetic reversals occur fully at random, with the distri- 
bution of jumping times following a Foisson distribution. 
Any deviation from such distribution, for large N and 
large negative E values, signals the presence of invariant 
curves preventing the motion to switch from one peak to 
the other one. 

Interestingly, and somewhat reminiscent of a critical 
FT, as shown in Fig.|2|the magnetic-reversal time-scale 
diverges as a power law of E at the critical energy Emt- 
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E /N 

tnt 



E ^/N 

obs 



Stat 



E — Ef. 



N 



(7) 



where 7(A^) oc N. 

Remarkably, such a dynamical behavior can be quaH- 
tatively - and almost quantitatively - explained only by 
statistical considerations. Before that, let us remark that 
^ can be inverted for E at any fixed TV, thus giving a r- 
depending characteristic energy Erev (t) . We will see the 
physical meaning of Erev in a while; here we just note 
that, since Q holds only for E > Etnt, it will also be 

Thermo- Statistics: restricting attention to the mean 
field Hamiltonian (0), a detailed statistical analysis us- 
ing techniques from lar ge d eviations theory leads to def- 
inite predictions (see [l^ for details) concerning the 
microcanonical (i.e., purely state-counting at fixed E) 
probability distribution Pstatirriy, E), the topological en- 
ergy threshold Etnt, and the thermo-statistical energy 
threshold Egtat, where the system undergoes a contin- 
uous ferro/paramagnetic PT in the microcanonical de- 
scription. 

Dynamics: Adopting a viewpoint as close as possible 
to the experimental one, we introduce an observational 
time tabs, during which the experiment can be performed. 
Assuming the experimentally measured value and the dy- 
namically observed time-average to coincide, we intro- 
duce the tofjs-averaged magnetization: 



{my) 



dt my(t). 



(8) 



The quantity {my) as a function of the energy E can be 
compared with my{E), obtained from the most likely val- 
ues PstatiEjiJiy) (the maxima of the distribution). This 
is shown in Fig. ^ where we have also indicated the 
thresholds Estat and Etnt as vertical lines. 

Of course, below Etnt and above Estat the two aver- 
ages coincides, at variance with the region between them. 
Also, the effective transition given by {my) occurs at 
some intermediate energy, different from both Etnt and 
Estat- This is not completely surprising, since both the 
topological energy Etnt and the thermo-statistical energy 
Estat refer to some limiting properties, the former when 
i — > cxD at any finite N, the latter when N oo at any 
time t. It is interesting to note that, inverting Eq. (TJ, to 
the observation time tabs there correspond a dynamical 
energy Erev{tobs) which is in fairly good agreement with 
the observed transition energy, see Fig. 0] 

Note also that, to a rather good accuracy 0,^3' there 
holds the following proportionality between the dynami- 
cal quantity Trev and the statistical quantity P{E, my): 



Trev ^ 



P 

Pq 




E/N 

FIG. 4: Numerical toba-averaged dynamical magnetization as 
a function of the specific energy e = E/N (circles). Dashed 
line are the analytical predictions from the thermo-statistical 
mean field model ©. Here is A*' = 60, a = 0, J? = 1. For the 
sake of comparison, both the thermo-statistical phase tran- 
sition and the topological non-connectivity thresholds have 
been indicated as vertical lines. Also indicated, as a vertical 
arrow, the expected dynamical magnetic-reversal threshold. 



where P = PstatiE, rhy) and Pq = Pstat{E, 0). Theoreti- 
cally ^3 , this can be traced back to the occurence of an 
entropy barrier AS* := S — Sq between the two compo- 
nents, where S = hiP and Sq — InPoi under the assump- 
tion of markovicity of the dynamical evolution. Detailed 
statistical consideration on the mean field model 10) lead 
to the analytical estimate jstat{N) — N to he compared 
with the numerical result jstatiN) « 0.88iV. 

Heuristically, ^ can be understood as follows. Sup- 
pose the dynamics is strongly enough chaotic, making the 
system almost freely wander through configuration space, 
in some sort of diffusive fashion. Then, the probability 
to encounter a configuration with my = 0, and therefore 
presumably switch from the my > to the the my < 
component, can be assumed to be simply proportional 
to the ratio of the number Mq of "right" configurations 
at m. 



y 



(9) 



to the number M+ of "wrong" configura- 
tions at my ^ 0, but counting only a half of the total 
wrong ones, namely those on the side where the system 
is coming from. (The system cannot know nothing yet 
about the other side still it's going to.) However, as long 
as the magnetization peaks are high and narrow, essen- 
tially all the "wrong" configuration are actually within 
the peak itself, i.e. normalizing by the total number 
M = M+ -I- Mq -\- M_ of energetically allowed config- 
urations gives the probabilities P+ ~ P. Finally, the av- 
erage time needed for magnetic reversal can be assumed 
to be taken proportional to the probability to encounter 
the right configuration with my = 0, times some char- 
acteristic "residence time" Tres during which the system 
essentially stays within each given iV-spin configuration. 
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along its longer wandering of the whole A'^-spin configu- 
ration space. Altogether, such a purely statistical argu- 
ment recovers exactly implicitly suggesting that the 
precise value of T^es might be the origin of the different 
coefficients found between the classical [Tsj and quantum 
|17| realizations of the same dynamical evolution equa- 
tions jnji- 

IV. CONCLUSIONS 

We showed the existence in classical Heisenberg spin 
models of a Topological Non-connectivity Threshold 
(TNT), caused by the anisotropy of the couphng when 
it induces an easy-axis of the magnetization. Below the 
TNT the constant energy surface is topologically discon- 
nected in two symmetrically equal components. 

Such a result on the Topology has connections both 



with the Dynamics (time-scales) and with the Thermo- 
Statistics (PT) of the system. In each component the 
magnetization along the easy axis has a definite sign, and 
it cannot change sign as time goes by, corresponding a fer- 
romagnetic behavior of the system. Above the TNT, in a 
strong chaotic regime, the magnetization reverses its sign 
with a characteristic time-scale which diverges as a power 
law at the TNT. Given enough chaos and enough obser- 
vational time, this leads to a paramagnetic behavior of 
the system. Moreover, the numerically observed link be- 
tween time-scales and probability distributions has been 
given an heuristical justification. 

The connection between the TNT and the range of 
the interaction has also been shown. For macroscopic 
spin systems the existence of this threshold determines a 
disconnected energy range that remains relevant (w.r.t. 
the total energy range) for long-range interactions, while 
it goes to zero for short-range interactions. 
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